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Abstract 

The latest CTEQ6.6 parton distributions, obtained by global analysis of hard scattering data in 
the framework of general-mass perturbative QCD, are employed to study theoretical predictions 
and their uncertainties for significant processes at the Fermilab Tevatron and CERN Large Hadron 
Collider. The previously observed increase in predicted cross sections for the standard-candle W 
and Z boson production processes in the general-mass scheme (compared to those in the zero-mass 
scheme) is further investigated and quantified. A novel method to constrain PDF uncertainties 
in LHC observables, by effectively exploiting PDF-induced correlations with benchmark standard 
model cross sections, is presented. Using this method, we show that the tt cross section can 
potentially serve as a standard candle observable for the LHC processes dominated by initial-state 
gluon scattering. Among other benefits, precise measurements of ti cross sections would reduce 
PDF uncertainties in predictions for single-top quark and Higgs boson production in the standard 
model and minimal supersymmetric standard model. 
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1. INTRODUCTION 



Parton distribution functions (PDFs) are essential inputs required to make theoretical 
predictions for the CERN Large Hadron Collider (LHC) and other hadron scattering facili- 
ties. They are extracted from a comprehensive global analysis of hard-scattering data from 
a variety of fixed-target and collider experiments in the framework of perturbative QCD. 
Experimental groups envision an ambitious program to tightly constrain PDF degrees of 
freedom using the upcoming LHC data. Such constraints will not be feasible in the early 
runs of the LHC, when experimental systematic errors will typically be large, and the collider 
luminosity itself will not be known to better than 10 — 20%. Thus the experiments plan to 
perform real-time monitoring of the collider luminosity through the measurement of bench- 
mark standard model cross sections, notably those for production of massive electroweak 
bosons [1-4]. These cross sections are large and can be measured fairly precisely soon after 
the LHC turn-on. To realize this goal, as well as to carry out the general physics program of 
the LHC, it is important, on one hand, to systematically explore the dependence of W. Z, 
and other "standard candle" cross sections on the PDFs and other aspects of QCD theory; 
and, on the other hand, to establish and exploit correlations with these observables arising 
from the dependence on the universal parton distributions. These tasks must be carried out 
using the most up-to-date PDFs, with quantitative estimates of their uncertainties. 

In a series of recent papers [5-7], we have extended the conventional CTEQ global PDF 
analysis [8, 9] to incorporate a comprehensive treatment of heavy-quark effects and to in- 
clude the most recent experimental data. The PDFs constructed in those studies consist 
of (i) the base set CTEQ6.5M, together with 40 eigenvector sets along 20 ortlionormal di- 
rections in the parton parameter space [5]; (ii) several PDF sets CTEQ6.5Sn (n=-2,...4), 
designed to probe the strangeness degrees of freedom under the assumption of symmetric or 
asymmetric strange sea [6]; and (iii) several sets CTEQ6.5XCn (n=0...6) for a study of the 
charm sector of the parton parameter space, in particular, the allowed range of independent 
nonperturbative ("intrinsic") charm partons in several possible models [7]. Some prominent 
physical consequences of these new PDFs compared to previous PDFs, particularly for W 
and Z production cross sections at hadron colliders, were pointed out in these papers. 

The current paper pursues a more detailed exploration of the physical implications of our 
new generation of PDFs, particularly at the LHC. In the process of this detailed study, we 
found it desirable to expand and improve the CTEQ6.5 analysis on several fronts. All new 
results presented in this paper are based on these improved PDFs, which we designate as 
CTEQ6.6.^ They will be described in Sec. 3. 

We make a systematic effort to address the quantitative challenges described in the first 
paragraph of this introduction. We study the impact of the new PDFs on the predicted 
cross sections for important physical processes at the LHC, with their associated uncer- 
tainty ranges. The PDF uncertainties in future measurements may behave at odds with 
initial intuitions because of rich connections between PDFs of different flavors and in dif- 
ferent kinematical ranges arising from (a) physical symmetries, such as scale invariance and 
parton sum rules, and (b) experimental constraints. In some cases, notably in Z boson 
production at the LHC, the largest PDF uncertainty arises from less-constrained subleading 



^ The CTEQ6.6 PDFs represent improvements over their counterparts in CTEQ6.0, while preserving the 
same physics inputs. Therefore we recommend that CTEQ6.6 be the preferred PDFs to use in future 
phenomenology studies. 
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scattering channels rather than from the well-known dominant subprocesses. In order to 
access such rich interconnections efficiently and completely, we introduce in Sec. 2 a quan- 
titative measure of correlations existing between the PDFs and physical observables. The 
correlation analysis extends the Hessian method [8, 10] to investigate pairwise relations be- 
tween collider observables. This analysis is employed in Sec. 4 as a tool to enhance the power 
of certain crucial phenomenological predictions, notably for W and Z boson production cross 
sections, and to examine the role of ti production as an additional standard candle process. 
At the end of Sec. 4, we apply the gained knowledge to identify the main sources of PDF 
uncertainty in single-top production and in several processes for Higgs boson production in 
the standard model (SM) and minimal supersymmetric standard model (MSSM). 

2. CORRELATIONS DUE TO PDFS 

In many applications, it is instructive to establish whether a collider observable shares 
common degrees of freedom with precisely measured SM cross sections through the nonper- 
turbative PDF parametrizations, a feature that can be exploited to predict the observable 
more reliably. In Section 4, we will explore such PDF-induced correlations between inter- 
esting collider cross sections. But first we will define the relevant theoretical framework. 

Let X be a variable that depends on the PDFs. It can be any one of the physical quantities 
of interest, or even a PDF itself at some given {x,ii). We consider X as a function of the 
parameters {aj} that define the PDFs at the initial scale /iq. Thus we have X{a), where a 
forms a vector in an iV-dimensional PDF parameter space, with being the number of free 
parameters in the global analysis that determines these PDFs. In the Hessian formalism for 
the uncertainty analysis developed in [10] and used in all of our recent work, this parton 
parameter space is spanned by a set of orthonormal eigenvectors obtained by a self-consistent 
iterative procedure [5, 8]. 

If Qq represents the best fit obtained with a given set of theoretical and experimental 
inputs, the variation of X{a) for parton parameters a in the neighborhood of Gq is given, 
within the Hessian approximation, by a linear formula 

AX{a) = X{a) - X{ao) = VXlg,, ■ Aa, (1) 

where VX is the gradient of X(a), and Aa = a — do. As explained in detail in Refs. [5, 8, 10], 
the uncertainty range of the PDFs in our global analysis is characterized by a tolerance factor 
T", equal to the radius of a hypersphere spanned by maximal allowed displacements Aa in the 
orthonormal PDF parameter representation. T is determined by the criterion that all PDFs 
within this tolerance hypersphere should be consistent with the input experimental data sets 
within roughly 90% c.l. The detailed discussions and the specific iterative procedure used 
to construct the eigenvectors can be found in Refs. [5, 8, 10]. 

In practice, the results of our uncertainty analysis are characterized by 2N sets of pub- 
lished eigenvector PDF sets along with the central fit. We have 2 PDF sets for each of the 
A^ eigenvectors, along the (±) directions respectively, at the distance |Aa| = T. The i-th 

— * 

component of the gradient vector VX may be approximated by 

f (2) 

where x\^^ and x\ ^ are the values of X computed from the two sets of PDFs along the (±) 
direction of the i-th eigenvector. The uncertainty of the quantity X due to its dependence 
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Figure 1: Dependence on the correlation ellipse formed in the 6X — 6Y plane on the value of cos (p. 
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where for simplicity we assume that the positive and negative errors on X are the same.^ 

We may extend the uncertainty analysis to define a correlation between the uncertainties 
of two variables, say X(a) and Y{a). We consider the projection of the tolerance hypersphere 
onto a circle of radius 1 in the plane of the gradients VX and VY in the parton parameter 
space [10, 11]. The circle maps onto an ellipse in the XY plane. This "tolerance ellipse" is 
described by Lissajous-style parametric equations. 



X = Xo + AXcos^, 

Y = Yo + AYcos{e + ^), 



(4) 
(5) 



where the parameter 9 varies between and 27r, Xq = X(ao), and Yq = Y{do). AX and AY 
are the maximal variations 6X = X — Xq and 6Y = Y — Yo evaluated according to Eq. (3) , 
and if is the angle between VX and VY in the {aj} space, with 



VX- vy 
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(6) 



The quantity cos (p characterizes whether the PDF degrees of freedom of X and Y are 
correlated {cos(p ~ 1), anti-correlated {cosip ~ —1), or uncorrelated (cosLp ^ 0). If units 
for X and Y are rescaled so that AX = AY (e.g., AX = AY = 1), the semimajor axis of 
the tolerance ellipse is directed at an angle 7r/4 (or 37r/4) with respect to the AX axis for 
cosip > (or cos(p < 0). In these units, the ellipse reduces to a line for cosip = ±1 and 



^ A more detailed equation for AX accounts for differences between the positive and negative errors [9, 11]. 
It is used for tt cross sections in Table 2 and Fig. 12. 
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becomes a circle for coscp = 0, as illustrated by Fig. 1. These properties can be found by 
diagonalizing the equation for the correlation ellipse, 



A magnitude of | cos 9?| close to unity suggests that a precise measurement of X (constrain- 
ing 5X to be along the dashed line in Fig. 1) is likely to constrain tangibly the uncertainty 
5Y in Y, as the value of Y shall lie within the needle-shaped error ellipse. Conversely, 
cos ^ implies that the measurement of X is not likely to constrain SY strongly.^ 

The parameters of the correlation ellipse are sufficient to deduce, under conventional 
approximations, a Gaussian probability distribution P(X, F|CTEQ6.6) for finding certain 
values of X and Y based on the pre-LHC data sets included in the CTEQ6.6 analysis. If 
the LHC measures X and Y nearly independently of the PDF model, a new confidence 
region for X and Y satisfying both the CTEQ6.6 and LHC constraints can be determined 
by combining the prior probability P{X, F|CTEQ6.6) with the new probability distribution 
P(X, y|LHC) provided by the LHC measurement. For this purpose, it suffices to construct 
a probability distribution 

P{X, F|CTEQ6.6+LHC) = P{X, F|CTEQ6.6)P(X, F|LHC), (8) 

which establishes the combined CTEQ6.6+LHC confidence region without repeating the 
global fit. 

The values of AX, AY, and cos(f are also sufficient to estimate the PDF uncertainty of 
any function f{X, Y) of X and Y by relating the gradient of f{X, Y) to dxf = df/dX and 
dyf = df/dY via the chain rule: 

A/ - |V/| - ^J {AX dxf f + 2 AX AY cosip dxf dyf + {AY dyff. (9) 

Of particular interest is the case of a rational function f{X,Y) = X'^/Y'^, pertinent to 
computations of various cross section ratios, cross section asymmetries, and statistical sig- 
nificance for finding signal events over background processes [11]. For rational functions 
Eq. (9) takes the form 



Af [AXy AX AY f \ 

For example, consider a simple ratio, / = X/Y. Then A///o is suppressed (A///o ~ 
|AX/Xo — AY/Yo\) if X and Y are strongly correlated, and it is enhanced (A///o ~ 
AX/Xq + AY/Yo) if X and Y are strongly anticorrelated. 

As would be true for any estimate provided by the Hessian method, the correlation angle 
is inherently approximate. Eq. (6) is derived under a number of simplifying assumptions, 
notably in the quadratic approximation for the function within the tolerance hypersphere, 
and by using a symmetric finite-difference formula (2) for {diX} that may fail if X is not 



^ The allowed range of SY/AY for a given 6 = SX/AX is ry ^ < SY/AY < where = 6cos(p ± 

Vl — S'^ sin</j. 
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monotonic. With these limitations in mind, we find the correlation angle to be a convenient 
measure of interdependence between quantities of diverse nature, such as physical cross 
sections and parton distributions themselves. For collider applications, the correlations 
between measured cross sections for crucial SM and beyond SM processes will be of primary 
interest, as we shall illustrate in Sec. 4. As a first example however, we shall present some 
representative results on correlations between the PDFs in the next section. 

3. OVERVIEW OF CTEQ6.6 PDFS 
A. CTEQ6.6 versus CTEQ6.1 

The CTEQ6.5 PDFs [5-7] and their improved version CTEQ6.6 presented here are based 
on a new implementation of heavy-quark mass effects in perturbative QCD cross sections, 
realized in the ACOT general-mass (GM) variable number scheme [12, 13] and supplemented 
by a unified treatment of both kinematical and dynamical effects according to the modern 
SACOT [12, 14] and ACOT-x [15] concepts. This improvement leads to significant changes 
in some key predictions compared to the ordinary zero-mass (ZM) scheme. The quality of 
the global analysis is further enhanced by to the inclusion of newer data sets, replacement 
of structure functions F2{x,Q) and F^,{x,Q) by the directly measured cross sections in HI 
and ZEUS deep-inelastic scattering (DIS) data sets, relaxation of ad hoc constraints on the 
parametrization of strange quark PDFs, and improvements in the global fitting procedure. 
The CTEQ6.6 PDF set includes four additional PDF eigenvectors to accommodate the 
free strangeness parametrization, as described below. It is also in a better agreement with 
HERA charm production cross sections than CTEQ6.5. The public distributions, available 
at the LHAPDF depository [16], include a central PDF set, denoted as CTEQ6.6M, and 44 
eigenvector sets that span the range of uncertainties in the parton parameter space due to 
input experimental errors. 

We illustrate the impact of the improved treatment of heavy-quark scattering by compar- 
ing predictions made using the CTEQ6.6 PDFs (realized in the GM scheme) and zero-mass 
CTEQ6.1 PDFs [9] (realized in the ZM scheme). All of the cross sections in our global anal- 
ysis are calculated at next-to-leading order (NLO) in perturbative QCD. Figures showing 
comparison of CTEQ6.6 and CTEQ6.1 for various PDF flavors are collected at [17]. Since 
these figures are rather similar to their counterparts comparing the CTEQ6.5 and CTEQ6.1 
PDFs [5], we do not reproduce them in this paper, except for the strangeness and charm 
PDFs. 

Incorporation of the full heavy-quark mass effects in the general-mass formalism leads 
to the suppression of heavy-flavor contributions to the DIS structure functions F\{x,Q) 
compared to the zero-mass formalism. For neutral-current Fx{x,Q), which dominate the 
global analysis, the suppressions occur in (a) the flavor-excitation partonic processes with 
incoming c and b quarks, through the rescaling of the light-cone momentum fraction variable; 
and (b) the light-flavor scattering processes involving explicit flavor creation (such as the 
gluon-photon fusion), through the mass dependence in the partonic cross section (Wilson 
coefficient) [5, 15]. Since the theoretical calculations in the global fit must agree with the 
extensive DIS data at low and moderate Q, the reduction in c, b, and g contributions in 
the GM NLO fit must be compensated by larger magnitudes of light quark and antiquark 
contributions. One therefore expects to see an increase in the light-quark PDFs extracted 
from CTEQ6.6 compared to those from CTEQ6.1 analyses in the appropriate {x, Q) region. 
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Figure 2: CTEQ6.6 PDF uncertainty bands (green shaded area) and CTEQ6.1M PDF (red solid 
line) for s = s at factorization scales /i = 3.16 and 100 GeV. 

The most consequential differences between CTEQ6.6 and CTEQ6.1 PDFs for u and d 
quarks occur at x < 10~^, cf. Ref. [17]. They can substantially affect predictions for quark- 
induced processes at the LHC, making the predicted cross sections larger by several percent 
(6-7% in W, Z production), as will be discussed in more detail in Sec. 4. 

As a new feature, the CTEQ6.6 analysis allows the shape of strange quark distributions 
to be independent from the non-strange sea distributions. We no longer impose the familiar 
ansatz s{x, /xq) oc u{x, fio) + d{x, fio), because the included dimuon DIS data {uA — > fi~^fi~X) 
[18] probes the strange quark distributions via the underlying process sW c, making 
the above ansatz unnecessary. However, as shown in Ref. [6], the existing experimental 
constraints on the strange PDFs remain relatively weak and have power to determine at most 
two new degrees of freedom associated with the strangeness in the limited range x > 10^^. 
Thus in the CTEQ6.6 analysis we add two new free parameters characterizing the strange 
PDFs.'' We continue to assume s{x) = s{x) in these fits, since the current data do not place 
statistically significant constraints on the difference between s{x) and s{x) [6]. 

At X < 10^^, the available data probes mostly a combination (4/9) [u{x) + u{x)] + 
(1/9) [d{x) + d{x) + s{x) + s{x)~\ accessible in neutral-current DIS, but not the detailed 
flavor composition of the quark sea. The shape of s{x,Ho) at very small x can vary 
over a large range, accompanied by corresponding adjustments in the other sea quark 
flavors. In other words, the strangeness to non-strangeness ratio at small x, Rg = 
lim^^o [^(x, /io)/ (u{x,fio) + J(x, yUo))], is entirely unconstrained by the data. But, on gen- 
eral physics grounds, one would expect this ratio to be of order 1 (or, arguably, a bit smaller). 



Assignment of more than two free strangeness parameters does not tangibly improve the quality of the fit 
and creates undesirable flat directions in the Hessian eigenvector space [6]. 
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Figure 3: CTEQ6.6 PDF uncertainty band for c = c with radiatively generated charm only (green 
shaded area) and PDFs with intrinsic charm: BHPS form with moderate (long-dashed) and strong 
(long-dash-dotted) IC; sea-like form with moderate (short-dashed) and strong (short-dash-dotted) 
IC. The red sohd line is for the CTEQ6.1M PDF. A factorization scale // = 100 GeV is assumed. 

Thus, in the current CTEQ6.6 analysis, we adopt a parametrization for the strange PDF of 
the form s{x,fio) = Aqx"^^ (1 — x)'^^P{x), where Ai is set equal to the analogous parameter 
of u and d based on Regge considerations. A smooth function P{x) (of a fixed form for 
all 45 CTEQ6.6 PDF sets) is chosen to ensure that the ratio Rs stays within a reasonable 
range. There is considerable freedom in the choice of P{x); and that is a part of the theo- 
retical uncertainty associated with any parametrization of the initial PDFs.^ Fig. 2 shows 
the CTEQ6.6 uncertainty bands of s{x,fi) at two values of /i, along with the corresponding 
CTEQ6.1M PDFs. The values of Rg for the 44 eigenvector sets in this PDF series span 0.63 
- 1.15. 

Since the PDFs are de facto used for ultra-high energy and astrophysics applications at 
X < 10~^, the CTEQ6.6 PDFs are tabulated down to x = 10~^ to provide numerically stable 
PDF values obtained by QCD evolution from parametrized initial parton distributions in 
this extreme x region. At the initial scale /xq, the PDFs are extrapolated into the region 
X < 10~^, not covered by the fitted data, by a Regge-like functional form proportional to 
x", where the negative parameter a is found from the fit at x > 10~^. At larger /i scales, 
the PDFs are predicted from the initial condition at /i = /io based on the NLO DGLAP 
evolution. No extra dynamical effects that may be significant at small x are included in this 
parametrization. 



^ The theoretical uncertainty associated with the small-x quark flavor composition is generally larger than 
the Hessian uncertainty bands due to propagation of experimental errors for a given P{x). 
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B. Fits with nonperturbative charm 

In the general-purpose CTEQ6.6 PDFs, we assume there is no nonperturbative intrinsic 
charm (IC), so that c(a;, /xq) = c{x,no) = at the initial evolution scale fiQ = rric — 1.3 GeV. 
To facilitate studies of the effect of possible IC, we have also created fits in which various 
amounts of IC (with c = c) are assumed. These fits amount to updated versions of the light- 
cone motivated (BHPS [19]) and the sea-like models for the shape of the input charm PDF 
discussed in [7]. For each IC model, we provide two PDF parameterizations with moderate 
and strong IC contributions, corresponding to charm quarks carrying 1% and 3.5% of the 
parent nucleon's momentum at /j, — rric, respectively. Figure 3 shows that the assumption of 
IC can substantially increase the amount of c = c at factorization scales as large as 100 GeV. 

C. Correlations 

As mentioned before, correlations can be computed for any variables, including the PDFs 
themselves. The x and /i dependence of the PDF-PDF correlations provides broad insights 
about theoretical and experimental constraints affecting physical cross sections. To explore 
this dependence, we present contour plots of the correlation cosine cosip for pairs of PDFs 
/oi {xi, Hi) and /a2(a^2, IJ'2) at scales = 2 and 85 GeV, plotted as a function of momentum 
fractions Xi and X2. A few such plots relevant for the ensuing discussion are shown in Fig. 4. 
A complete set of the contour plots in color and grayscale versions is available at [20]. 

Light (dark) shades of the gray color in Fig. 4 correspond to magnitudes of cos (/? close to 1 
(-1), as indicated by the legend. Several interesting patterns of the PDF correlations can be 
generally observed. First, consider self-correlations, in which the correlation cosine is formed 
between two values of the same PDF (oi — a2, iJ>i — evaluated at momentum fractions 
xi and X2. The examples include the u — u,g — g, and c — c correlations at iJ>i — 112 — 85 GeV 
shown in Figs. 4(a)-(c). 

Each self-correlation plot includes a trivial correlation, cos 93 ~ 1, occurring when Xi and 
X2 are about the same. This correlation occurs in light-colored areas along the xi — X2 
diagonals in Figs. 4(a)-(c), of the shape that depends on the flavor of the PDF and the 
associated scale. 

In the case of an up quark [Fig. 4(a)], the trivial correlation is the only pronounced pattern 
visible in the contour plot. The gluon PDF [Fig. 4(b)] and related c, b PDFs [Fig. 4(c)] show 
an additional strong anti-correlation in the vicinity of {xi,X2) ~ (0.2,0.01) arising as a 
consequence of the momentum sum rule. Important implications of this anticorrelation will 
be discussed in Section 4. 

The gluon-strangeness correlation in Fig. 4(d) illustrates some typical patterns encoun- 
tered in the case of sea partons. The gluons show a correlated behavior with strange quarks 
(and generally, sea quarks) at a; < 10~^ — 10~^ as a reflection of the singlet evolution (a light 
area in the lower left corner) . The gluon PDF is anticorrelated with the strangeness PDF at 
Xi ~ 0.1, a;2 < 0.01 because of the momentum sum rule (a dark area at the bottom). This 
anticorrelation significantly affects predictions for W and Z cross sections at the LHC. More 
complicated (and not so well-understood) patterns occur at 2 > 0.3, where the sea-parton 
behavior is less constrained by the data. Other types of correlation patterns, associated with 
the sum rules, perturbative evolution, and constraints from the experimental data, can be 
observed in the full set of contour plots [20]. 
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Figure 4: Contour plots of correlations cost^ between two PDFs and /2(x2,/i) at /x = 

85 GeV: (a) u — u; {h) g — g; (c) c — c; (d) s — g. Both axes are scaled as x^'^. 

4. IMPLICATIONS FOR COLLIDER EXPERIMENTS 

Soon after its turn-on, the LHC will provide vast samples of data in well-understood 
scattering processes at the electroweak scale, notably production of massive weak bosons 
and Z^. These data will facilitate useful experimental calibrations as well as measurements 
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of the LHC luminosity and PDFs with tentative experimental accuracy of about 1% [1-3]. 
To utilize such measurements of "standard candle" cross sections most productively, one 
must understand how they constrain PDF degrees of freedom. It is necessary to explore 
how the predicted cross sections change due to improvements in the theoretical model, such 
as the transition from the zero-mass to general-mass factorization scheme, as well as due to 
the remaining freedom in the PDF parameters allowed by the data in the global fit. 

In this section, we examine the quantitative connections between the PDFs and physical 
observables by applying the correlation analysis introduced in Section 2. Since only large 
correlations (or anticorrelations) will be useful for practical purposes, we concentrate on 
pairs of cross sections characterized by large magnitudes of the correlation cosine, chosen in 
this paper to be |cos(/?| > 0.7. Overall dependence on the heavy-quark scheme is explored 
in Section 4 A. A focused study of total inclusive cross sections is presented for W, Z, top- 
quark, and Higgs boson production in Sections 4B, 4C, and 4D. The total cross sections 
are computed at the next-to-leading order (NLO) in the QCD coupling strength ctg, using 
W boson mass Mw — 80.403 GeV, Z boson mass Mz — 91.1876 GeV, and top quark mass 
mt = 171 GeV. 

In all calculations, both the renormalization and factorization scales are set to be equal to 
the mass of the final-state heavy particle, unless specified otherwise. The results presented 
below are representative of generic patterns observed in the PDF dependence of the studied 
cross sections. They depend weakly on theoretical assumptions about the quark fiavor 
composition at small x discussed in Section 3. 

The LHC collaborations intend to measure ratios of LHC cross sections to "standard- 
candle" (especially W and Z) cross sections. If two cross sections X and Y share common 
systematics, both experimental and theoretical, the systematic imcertainties partially cancel 
in the ratio r — X/Y . This cancellation is especially important in the first year or two of 
the LHC running, when the collider luminosity will only be known to the order of 10 — 20%. 
Eventually the LHC ratios can provide important precision tests of the Standard Model, 
such as a precise measurement of the W boson mass from the ratio of W and Z total cross 
sections [21-24]. 

The PDF uncertainty on r is reduced (Ar/r ~ \AX/X — AY/Y\) if X and Y are strongly 
correlated, cf. Eq. (10). It is enhanced {Ar/r ~ AX/X + AY/Y) ii X and Y are strongly 
anticorrelated. It is therefore beneficial to construct the cross section ratios from pairs of 
correlated cross sections to reduce the PDF uncertainties. We will examine a ratio between 
the highly-correlated W and Z total cross sections in Section 4B3. We will also explore 
the inclusive tt production cross section, potentially useful for building ratios when an anti- 
correlation with W and Z cross sections occurs. W, Z, tt total cross section values for 45 
CTEQ6.6 sets suitable for the computation of various PDF correlations can be downloaded 
at [17]. 

Leading channels for neutral Higgs boson (/i°) production in the standard model [25] are 
also investigated, including gluon-gluon fusion gg — > the dominant Higgs production 
mechanism; vector boson fusion qq WW — > /i° (or qq — > ZZ h^), a prominent 
channel employed at the LHC for both discovery and measurement of h^WW (or hPZZ) 
couplings; and associated production of Higgs boson with a massive weak boson, qq —>■ Wh^ 
{qq — > the leading discovery mode for Higgs masses nih < 200 GeV at the Tevatron. 

In addition, Higgs production channels in MSSM are considered: neutral CP-even or 
CP-odd Higgs boson production via bottom-quark annihilation, bb ^ h {h = h^, H^, 
[26]; charged Higgs boson production via s, c, b scattering, cs + cb — > h'^ [27]; and associated 
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production qq' W ^ Ak^ of CP-odd {A) and charged Higgs bosons [28, 29], with Higgs 
masses related as = m\ + at the Born level. The bb ^ h and cs + cb ^ cross 
sections may be greatly enhanced if the ratio tan/3 = Vu/va of vacuum expectation values 
for up- and down-type Higgs doublets is of order 10 or more, a possibility that remains 
compatible with LEP, Tevatron, and other constraints [30]. To be specific, we evaluate qqh 
couplings at the tree level and choose effective tan/? = 50; but most of our results (presented 
as cross section ratios for different PDF sets and correlation cosines) are independent of the 
tan (3 value. 

Supersymmetric neutral Higgs production 66 — > /i is sensitive largely to b{x, fj,) and 
g{x, n), while charged Higgs production probes a combination of s{x, fj,), c{x, fj,), and b{x, /i). 
Long-distance hadronic functions entering these processes require a general-mass scheme 
approach [6, 7, 31, 32]. In hard-scattering inclusive cross sections, heavy-flavor mass correc- 
tions enter through ratios mlj^/pf and are suppressed at large momentum scales ^ nT'l j^. 
We therefore neglect the charm and bottom quark masses in hard-scattering matrix elements 
at TeV energies, while keeping the mass effects inside the PDFs. 

A. Total cross sections at the Tevatron and the LHC 

To discuss common differences between the general-mass (CTEQ6.6) and zero-mass 
(CTEQ6.1) predictions, we compute several total cross sections sensitive to light-quark, 
gluon, or heavy-quark scattering using the NLO programs WTTOT [33] and MCFM [34- 
36]. Figure 5 shows several CTEQ6.6 and CTEQ6.1 total cross sections and their PDF 
uncertainties at the pp collider Tevatron {^/s = 1.96 TeV) and the pp Large Hadron Collider 
(a/s — 14 TeV). They are plotted as ratios to CTEQ6.6M predictions. We do not show 
predictions for the CTEQ6.5 set, as those agree with CTEQ6.6 within the uncertainties. 

The processes shown in the figure for the Tevatron occur at relatively large momentum 
fractions x, where CTEQ6.6 and CTEQ6.1 PDFs agree well. Consequently the CTEQ6.6 and 
CTEQ6.1 cross sections for the Tevatron processes coincide within the PDF uncertainties, 
as illustrated in Fig. 5(a). The magnitudes of the PDF uncertainties also remain about the 
same. 

At the LHC [Fig. 5(b)], CTEQ6.6 cross sections for W and Z boson production are 
enhanced by 6 — 7% compared to CTEQ6.1 because of the larger magnitudes of u \x,Q) 
and (i:\x,Q) in the relevant range of a; ~ Mw,z/V^ = 10"^ - 10"^. The CTEQ6.6 light- 
quark parton luminosities Cq^q.{xl,X2, Q) = qi{xi, Q) qj{x2, Q) (and therefore CTEQ6.6 cross 
sections) are larger at such x hy 2Sf ~ 6%, where Sf ~ 3% is the typical increase in 
the CM light-quark PDFs compared to the ZM PDFs. The Hessian PDF error obtained 
by our standard 90% c.l. criterion has decreased from 4.5 — 5% in CTEQ6.1 to 3.5% in 
CTEQ6.6, mostly because new DIS experimental data were included in the CTEQ6.6 fit. The 
differences between CTEQ6.6 and CTEQ6.1 exceed the magnitude of the NNLO correction 
to W and Z cross sections of order 2% [37-40], indicated by a dashed line. Uncertainties of 
this size have important implications for the calibration of the LHC luminosity. 

Other cross sections dominated by light (anti)quark scattering at a; ~ 10^^ increase by 
comparable amounts. For example, the CTEQ6.6 cross sections for associated W^h or 
Z^h boson production exceed CTEQ6.1 cross sections by 3 — 4%. In contrast, in processes 
dominated by gluon or heavy-quark scattering, such as fl production or gg h^X, the 
general tendency for the CTEQ6.6 cross sections is to be a few percent lower compared 
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Figure 5: Representative CTEQ6.6 (red boxes) and 6.1 (blue stars) total cross sections and PDF 
uncertainties at the Tevatron and LHC, normalized to the CTEQ6.6M cross section. Green triangles 
indicate CTEQ6.6 cross sections obtained under an assumption of a strong sea-like intrinsic charm 
production (IC-Sea-3.5%). 



to CTEQ6.1. The CTEQ6.6M cs + cb cross section is enhanced with respect to 

CTEQ6.1M by its larger strangeness PDF (cf. Fig. 2), despite some suppression of the cb 
contribution to this process. The LHC cross sections in Fig. 5(b) may change substantially if 
a fraction of charm quarks is produced through the nonperturbative "intrinsic" mechanism. 
For example, the rate for production sc + sb ^ of MSSM charged Higgs bosons with 
mass 200 GeV would increase by 30% if sea-like intrinsic charm contributions carry 3.5% of 
the parent nucleon's momentum, the maximal amount tolerated in the fit. 
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W* & Z cross sections at the Tevatron Run-2 
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t igure 6: CTEQ6.6 NNLL-NLO TF± and Z production cross sections in the Tevatron Run-2, 
compared with predictions for other PDF sets by (a) CTEQ and (b) other groups. 



B. W and Z boson production cross sections 

1. CTEQ6.6 vs. other PDF sets 

Large groups of collider cross sections, notably those dominated largely by quark scat- 
tering or largely by gluon scattering, exhibit correlated dependence on PDFs. Production 
of charged and neutral bosons are essential quark-quark scattering processes that 
show such correlations. 

Figs. 6 and 7 summarize predictions for the total cross sections a ofW and Z production, 
obtained in the NNLL-NLO resummation calculation (of order C(Q;s)-|-leading higher-order 
logarithms) [41-43] and using recent PDF sets by CTEQ [5, 8, 9, 44-48], Alekhin/AMP 
[49, 50], HI [51], MRST/MSTW [52-55], and ZEUS [56, 57] groups. We include decays of 
the massive bosons into lepton pairs in the improved Born approximation [42, 58]. To provide 
a visual measure of the PDF uncertainty, each figure shows an error ellipse corresponding 
to our usual tolerance criterion.^ The ellipses are found from Eqs. (4) and (5), using the 
parameters listed in Table 1. 

At the Tevatron, the CTEQ6.6 cross sections lie close to the CTEQ6.1 and CTEQ6.5 
predictions [Fig. 6(a)], i.e., the dependence on the heavy-fiavor scheme is relatively weak in 
this case. The CTEQ6.6 cross sections agree well with predictions based on the PDF sets 
by the other groups [Fig. 6(b)]. 

At the LHC, the general-mass CTEQ6.6 cross sections exceed the zero-mass CTEQ6.1 
cross sections by 6 — 7%, as shown in Fig. 7(a). The CTEQ6.6 and 6.5 Z^ cross sections are 



^ In two-dimensional plots, this criterion corresponds to probability somewhat smaller than 90%. 
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r igure 7: CTEQ6.6 NNLL-NLO and Z production cross sections in the LHC, compared with 
predictions for (a) different PDF sets by CTEQ and (b) other groups. 



about the same. The CTEQ6.6 cross section is somewhat smaller than the CTEQ6.5 
cross section and does not lie on the same line in the W — Z plane as the previous CTEQ 
sets. The predictions based on the latest CTEQ6.6, MSTW'06, and AMP'06 PDFs agree 
within 3% [Fig. 7(b)]. 

The total cross sections shown here are somewhat affected by higher-order contributions, 
not included under the current NNLL-NLO approximation. A prediction of absolute magni- 
tudes of W and Z cross sections with accuracy 1% would require to simultaneously evaluate 
NNLO QCD contributions of order 0{a1) [37-40] and NLO electroweak contributions of or- 
der 0{asCtEw) [59-67] for both the hard cross sections and PDFs, including all relevant spin 
correlations [68, 69]. This level of accuracy is not yet achieved. However, the higher-order 
terms rescale the NLO hard cross sections by overall factors with weak dependence on the 
PDFs [39, 40]. Therefore our NLO-NNLL total cross sections should reasonably estimate 
true relative differences caused by the PDFs. 

The NNLL-NLO cross sections for AMP'2006 and MSTW'2006 are computed using the 
NNLO PDFs, since the NLO PDFs for these sets are not available. Such combination is 
acceptable at the order we are working. For example, mixing of NLO-NNLL and NNLO 
orders is of little consequence in the case of the MRST'2004 set, for which the replacement of 
the NLO PDFs by the NNLO PDFs changes the cross sections by about 1.5%. Variations of 
this magnitude are clearly permissible until a more precise computation is fully developed. 



2. Correlations between W, Z cross sections and PDFs 

Although strong PDF-induced correlations between the W and Z cross sections are ob- 
served at both colliders, the mechanism driving these correlations at the LHC is not the 
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same as at the Tevatron. The essential point is that, although the weak bosons are mostly 
produced in u and d quark-antiquark scattering, this dominant process may contribute little 
to the PDF uncertainty because of tight constraints imposed on the up- and down-quark 
PDFs by the DIS and Drell-Yan data. 

Instead, a substantial fraction of the PDF uncertainty at the LHC (but not at the Teva- 
tron) is contributed by sizable, yet less constrained, contributions from heavy-quark (s, c, 
h) scattering. Subprocesses with initial-state s, c, and, to a smaller extent, h and g partons 
deliver up to 20% of the NLO rate at the LHC, compared to 2-4% at the Tevatron. All 
these partons are correlated with the gluons via DGLAP evolution, so that the LHC W 
and Z cross sections are particularly sensitive to the uncertainty in the gluon PDF. Conse- 
quently, W and Z cross sections at the LHC are better correlated with processes dominated 
by gluon scattering in the comparable kinematical range and not necessarily with u and d 
quark scattering, in striking contrast to the Tevatron. 

To illustrate this point. Fig. 8 shows correlation cosines (cos<^) between the W, Z cross 
sections and PDFs fa{x,Q) of different flavors, evaluated as functions of the momentum 
fraction a; at Q = 85 GeV. The largest correlations between the cross section and PDFs occur 
at momentum fractions x of order M\y,z/\/s corresponding to central rapidity production, 
i.e., at X ~ 0.04 at the Tevatron and 0.006 at the LHC. PDF flavors with a very large 
correlation are associated with the major part of the PDF uncertainty in the physical cross 
section. Additional constraints on this flavor would help reduce the PDF uncertainty. 

At the Tevatron Run-2 [Figs. 8(a) and 8(b)], large correlations exist with u, u, d, and 
d PDFs, with cosip reaching 0.95. No tangible correlation occurs with PDFs for s, c, b 
(anti)quarks and gluons. 

At the LHC [Figs. 8(c) and 8(d)], the largest correlations are driven by charm, bottom, 
and gluon PDFs, followed by smaller correlations with u, d, and s quarks. The correlation 
with the u and d PDFs is reduced, although not entirely eliminated. A large correlation with 
the gluon at a; ~ 0.005 is accompanied by a large anti-correlation (cost/? ~ —0.8) with the 
gluon at a; ~ 0.1 — 0.2, as a consequence of nucleon momentum conservation (cf. Section 3C). 
This feature implies that the LHC W. Z cross sections are strongly anticorrelated with 
new particle production in gluon or heavy-quark scattering in the TeV mass range, and 
moderately anticorrelated with tt production (cf. Section 4C2). 

3. Impact of strangeness 

We stress that the large correlation between the LHC Z cross section and the gluon 
distribution is generated primarily by the sizable ss scattering contribution, and not by a 
smaller qg scattering contribution. Similarly, a large correlation in W boson production 
is driven for the most part by the cs contribution. Yet the correlations of Z and W cross 
sections with the strangeness PDF, taken on its own, are rather small [see Figs. 8(c) and (d)]. 
Adjustments in s{x, Q) can be easily compensated by other quark PDFs without changing 
the physical W and Z cross sections. This extra freedom is absent in the case of the gluon 
PDF, since it affects all sea PDFs at once through perturbative evolution. 

On the other hand, the ratio rzw = (^z/icw+ +c"vf-) of the Z° and cross sections is 
very sensitive to the uncertainty in strangeness. Nominally rzw is an exemplary "standard 
candle" LHC observable because of the cancellation of uncertainties inside the ratio. The 
CTEQ6.6 prediction rzw = 0.100 ± 0.001 is in an excellent agreement with the predictions 
based on the other recent PDF sets. This result can be derived from Eq. (10) by substituting 
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Figure 8: (a,b) Correlation cosine (cosc/j) between the total cross sections for and production 
at the Tevatron and PDFs, plotted as a function of a; for Q = 85 GeV; (c,d) the same for the LHC. 



the correlated parameters of the W and Z cross sections in Table 1. 

Fig. 9 indicates that rzw is mostly correlated with the strangeness PDF s{x, Q) in the 
region 0.01 < a: < 0.05. It is anticorrelated with u ' and cT^ at x ~ 10~^. There is no tangible 
correlation with the gluon, charm, and bottom PDFs. 

Since the strangeness is the least constrained distribution among the light-quark fla- 
vors [6, 53, 70], it is thus important to correctly model its uncertainty to estimate rzw- 
CTEQ6.6 is our flrst general-purpose PDF set that includes an independent parametriza- 
tion for strangeness. It predicts a larger uncertainty in s{x, Q) (hence, a larger Ar^vi/) than 
the previous PDF analyses, which artiflcially linked s{x, Q) to the combination {u + S) of 
the lightest sea quarks. 

The increased fractional uncertainty ^rzwl^zw is related to weaker correlation between 
the CTEQ6.6 W and Z cross sections at the LHC. According to Eq. (10), ^rzw/^zw scales 
approximately as (1 — cosv^)^/^, given that the W and Z fractional uncertainties are about 
the same {^aw /o'w ~ ^^z/cy%)- The value of cosip decreases from 0.998 in CTEQ6.1 to 
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Figure 9: Correlation cosine (cos 99) between the ratio ozlf^w of LHC total cross sections for Z° 
and production at PDFs of various flavors, plotted as a function of x for Q = 85 GeV. 
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Figure 10: W and Z correlation ellipses at the LHC obtained in the fits with free and fixed 
strangeness, as well as with maximal intrinsic charm contribution. 

0.956 in CTEQ6.6, or 0.982 if s(x, Q) is fixed during the Hessian analysis at its best-fit 
CTEQ6.6M shape. As a result of the smaller cosy?, ISr zw I "f zw increases threefold in the 
CTEQ6.6 prediction, even though the fractional uncertainties A(Tvi/,z/c"vi/,z are reduced. 

The eccentricity of the oz — o'w tolerance ellipse grows with cos(p, implying a narrow 
CTEQ6.1 ellipse and a broader CTEQ6.6 ellipse shown in Fig. 10. The CTEQ6.6 ellipse 
narrows if the strangeness parameters are fixed at their best-fit values as described above. 
Very different values of rzw can be obtained if one allows for the "intrinsic charm" contribu- 
tion. In Fig. 10, the empty triangles correspond to the "maximal" intrinsic charm scenarios. 
These cross sections lie on the boundary or outside of the CTEQ6.6 tolerance ellipse. They 
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Figure 11: (a) CTEQ6.6 W~ and production cross sections at the LHC compared with pre- 
dictions for other PDF sets; (b) correlation (cos(/j) between the ratio a-sji/+ j Oyi^- of and 'W~ 
total cross sections and valence quark PDFs. 



can be potentially ruled out by measuring rzw precisely at the Tevatron, as can be deduced 
from Fig. 6(b). 



^. W vs cross sections 

The charge dependence of weak boson production at the LHC is explored by plotting 
W~ vs. cross sections in Fig. 11(a). In this case, the CTEQ6.6 prediction agrees v^ell 
v\rith the latest MSTW and AMP sets, while more substantial differences exist with the 
earlier MRST and ZEUS PDF sets. Possible intrinsic charm contributions (indicated by the 
IC-BHPS and IC-Sea points) would increase both cross sections. 

The cross section ratio aw-/o'w+ is most correlated with the valence w-quark PDF uy = 
u — u at Q = 85 GeV, followed by the valence d-quark [cf. Fig. 11(b)]. There is a large 
correlation with uv{x,Q) at a; ~ 0.1 and anticorrelation at x ~ 0.05. Other PDF flavors do 
not demonstrate pronounced correlation with aw-/crw+ and are not shown in the figure. 



C. Top-quark production and gluon uncertainty 

The prominent role of gluons in driving the PDF uncertainty in many processes has been 
pointed out in the past by noticing that most of the PDF uncertainty is often generated by a 
certain PDF eigenvector sensitive to the gluon parameters. While the abstract eigenvectors 
give largely a qualitative insight, the correlation analysis relates the PDF uncertainty directly 
to parton distributions for physical flavors at known {x, (i). We will now apply the correlation 
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technique to investigate the uncertainties associated with gluons and heavy quarks in two 
other prominent processes, production of top quark-antiquark pairs and single top quarks. 

1. Parametrization oftt and single-t total cross sections 

The inclusive rate for production of top quark-antiquark pairs, pp ' — > ttX, is measured 
with good precision in the Tevatron Run-2, and an even more precise measurement is possible 
with the large event yield expected at the LHC. Such a measurement can test the PQCD 
prediction and provide an alternative method to determine the mass of the top quark [71, 72]. 
At leading order, tt pairs are produced via qq scattering (contributing 85% of the rate at the 
Tevatron) and gg scattering (contributing 90% of the rate at the LHC). At NLO, both the 
fixed-order [73-76] and resummed [72, 77-88] cross sections have been computed. In this 
study, we calculate the NLO ti cross section using CTEQ6.6 PDFs and the MCFM code 
[34-36], for three values of the factorization scale (/U = 0.5 mt, nit, and 2mt). 

Production of single top quarks pp ' f^X [89-109] provides a unique means to measure 
the Wtb coupling with the goal of constraining new physics [91, 110-113]. It was recently 
observed for the first time by the Tevatron D0 [114] and CDF Collaborations [115]. We will 
focus on two single-top production channels, proceeding through t-channel and s-channel 
exchanges of charged weak (W) bosons. The t-channel W exchange involves bottom-quark 
scattering qb q't and dominates both the Tevatron and LHC rates. The s-channel W 
exchange is similar to conventional W boson production, but occurs at larger typical x 
values (of order mtj \/s rather than M^/y^). It may be observable at the Tevatron, but 
has a relatively small event rate at the LHC. We compute the NLO cross sections for single- 
top production using the programs from Refs. [103, 104] and [105-107]. 

For each scale we parametrize the resulting cross sections in the vicinity of the world- 
average top mass mt — 171 ± \.\{stai) ± 1.5(sysi) GeV [116] by a function 

a{mt, fi) = A{fi) + B{fi) {mt - 171) + C(/i) (m* - 171)' , (11) 

where the units of a and mt are picobarn and GeV. The variation with respect to the 
reference cross section (corresponding to /i = mt) gives the NLO scale dependence A^(mt), 
discussed in detail in Section 4D. The uncertainty in a due to the variation of within 
the experimentally allowed range gives the mass dependence, denoted by A^- We also 
calculate the relative PDF uncertainties, denoted by Ap£)p{mt). The values of A^, Ap^p, 
and coefficients A, B, C are listed in Table 2. In single-top cross sections, we set C = 0. A 
plot of this parametric dependence in tt production at the Tevatron and LHC is shown in 
Fig. 12. The correlation cosines between tt, single-top, W, and Z cross sections are listed in 
Tables 1 and 3. 

2. PDF-induced correlations 

The PDF dependence of the top-quark cross sections follows a few non-trivial trends, 
which can be understood by studying a;-dependent correlations between the top-quark cross 
sections and PDFs presented in Fig. 13. Our discussion will also refer to Tables 1-3. 

1. Contrary to the naive expectation, the main PDF uncertainty in tt production at 
the Tevatron is not associated with the leading qq scattering channel. Rather, the 
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mt[GcV] [GeV] 

(a) (b) 

Figure 12: CTEQ6.6 predictions for inclusive ti production at (a) the Tevatron Run-2 and (b) 
the LHC, showing the NLO cross section in pb versus the top-quark mass mt. The three curves 
correspond to three choices for the factorization scale: = mt/2, rrit, and 2mt. The error bars are 
the PDF uncertainties. Also shown is the la error range from the 2007 world-average experimental 
mt value by the Tevatron Electroweak Working Group, rrit = 171 ± 1.1 it 1.5 GeV [116]. 



uncertainty is mostly correlated with the gluon PDF probed at x ~ 0.3, as clearly 
shovi^n by the tt-PDF correlations in Fig. 13(a). At such x values, the quark PDFs are 
tightly constrained, resulting in a small uncertainty in the leading qq channel, while 
the gluons are poorly known, resulting in a very large uncertainty in the subleading 
gg channel. The net result is a substantial total PDF uncertainty, Apor = 7.4%, in 
the Tevatron ti cross section, contributed mostly by gg scattering. By the momentum 
sum rule, the Tevatron ti cross section is strongly anti-correlated with gluon scattering 
at a; ~ 0.05. It does not exhibit a strong (anti-)correlation with W or Z production at 
either collider, cf. Table 1. 

2. The Tevatron t-channel single-top cross section is mostly correlated with the 6-quark 
PDF at X ~ 0.2, as illustrated by Fig. 13(c). It has a substantial PDF uncertainty, 
ApDF =10.3%. It is correlated with the Tevatron ti cross section (cosy) = 0.81) and 
anticorrelated with Z, W production at the LHC {cosip = —0.82 and —0.79) through 
the shared correlation with the gluon at large x. 

3. At the LHC [Figs. 13(b) and (d)], the ti and t-channel single-top cross sections are 
also mostly correlated with g, c, and b PDFs, which, however, are well-constrained 
at X ~ 0.05 and 0.01, typical values in this case. The PDF uncertainty is of order 
3% in both processes. The LHC ti cross sections are anticorrelated with the Z and 
W cross sections. The ti — Z and ti — W correlation cosines are large and negative: 
COSV9 = —0.8 and —0.74, respectively. The strong anticorrelation reveals itself in the 
shape of the ti — Z ellipse, plotted in Fig. 14(a) by using the PDF error parameters 
in Table 1. A similar anticorrelation exists between the ti and W cross sections [117]. 
Interestingly enough, the LHC t— channel single-top cross section is only mildly anti- 
correlated with ti production and mildly correlated with W and Z production [see 
Table 3], in contrast to the Tevatron. 

4. At both colliders, the s-channel PDF uncertainty is of order 3% and correlated mostly 
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CTEQ6.6: correlation between o-,o, and f(x,Q-85. GeV) 
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Figure 13: PDF-induced correlations {cos ip) between the total cross sections for tt and single-top 
production, and PDFs of various flavors, plotted as a function of a; for Q = 85 GeV: Tevatron Run-2 
(left column); LHC (right column). 
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Z° vs tt production at the LHC at NLO 
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Figure 14: Correlation ellipse for tt and Z boson NLO total cross sections at the LHC. 

with M ' and cf' PDFs [cf. Figs. 13(e) and (f)], as well as with the Tevatron Z, W 
rates. Remarkably, the LHC s-channel single-top cross section is not particularly 
correlated with the LHC W, Z cross sections, despite its similarities with W boson 
production. This peculiarity is due to large x values (x ~ mt/y/s ~ 0.01) typical 
for the s-channel single-top production. At such x, charm and bottom initial-state 
contributions are relatively small and do not affect single-top production as much 
as W boson production, hence preventing the gluon-driven PDF uncertainty from 
contributing sizably to the s-channel cross section. 

5. The improved evaluation of heavy-quark terms in the CTEQ6.6 PDFs reduces the 
Tevatron (LHC) tt cross sections by 4% (3.5%) compared to CTEQ6.1. The CTEQ6.6 
Tevatron t-channel single-top cross section is about 6% smaller than the CTEQ6.1 cross 
section. The other three CTEQ6.6 single-top cross sections differ from the CTEQ6.1 
cross sections by less than 2%. 



m,=^=171 GeV 



CTEQ6.6^ 



D. tt production as a standard candle; Higgs boson production 

To recap the previous section, at the Tevatron tt production is strongly correlated with 
single-top production. At the LHC, it is strongly anticorrelated with Z, W production. 
Precise measurements of tt rates could provide valuable constraints on the gluon and heavy- 
quark PDFs, if the associated theoretical and experimental uncertainties are each reduced 
below 3 — 5%. These measurements would help bring down the PDF error in many processes 
and supply an alternative way to monitor the collider luminosity. If a cross section is anti- 
correlated with W and Z production, it could be normalized to the ti cross section. The 
PDF error will be suppressed in such a cross section ratio, in contrast to the ratio with the 
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W or Z cross section. Other systematic errors may cancel better, too, if the process shares 
common elements with ti production. 

Let us consider some specific examples. The PDF uncertainties for Z, W^, tt, and Higgs 
boson production via gluon fusion at the LHC, with Higgs mass rrih — 500 GeV, are 
3.4%, 3.5%, 3.2%, and 4%, respectively. The correlation cosines are 0.956 for Z and 
cross sections, 0.98 for and ti cross sections, and —0.87 for h'^ and Z cross sections. 
By Eq. (10) the PDF uncertainty on the Z/W^, h^/ti, /Z cross section ratios are 1.3%, 
1.5%, and 7.2%, i.e., the correlated cross sections produce ratios with the smallest PDF 
uncertainties. 

The viability of precise ti production measurements can be examined by studying the- 
oretical uncertainties on the ti cross section presented in Table 2. At present, the scale 
dependence is larger than the PDF uncertainty /S.pdf and top-mass uncertainty at 
both colliders, suggesting that higher-order (NNLO) corrections have a tangible impact on 
the ti rate.'' The NNLO contributions will be computed in the near future [118], which will 
likely reduce the associated uncertainty to a few percent. 

The top-mass uncertainty A^ can be brought down to 2-3% by measuring mt with ac- 
curacy of order 1 GeV, as planned by the Tevatron experiments. Further advancements can 
possibly improve accuracy in the treatment of heavy-quark mass effects at NNLO to about 
a percent level. The current (not related to the luminosity) experimental systematic error 
for the fi cross section at the Tevatron is 8%) [119], with further improvements likely. A 
similar systematic error of order 5% may be feasible at the LHC. 

The complementarity of constraints on Higgs boson searches from ti production at the 
LHC and Z boson production at the Tevatron and LHC is illustrated by Fig. 15, showing the 
correlation cosines (cos </?) between Higgs, Z, and ti production cross sections as a function of 
Higgs mass. These cosines are shown by lines, in addition to separate markers corresponding 
to correlations between W, Z, and top production processes discussed above. The relevant 
cross sections are collected in Table 1 for Z, and ti production, and in Table 4 for Higgs 
production processes. 

As discussed earlier, there is a very large correlation between Z boson, boson, and 
boson production at the LHC, and a strong anticorrelation between Z boson and ti 
production. Only a mild correlation exists between Z production at the Tevatron and Z 
production at the LHC. There is a moderate correlation between Z production and the 
production of a light (120 GeV) Higgs boson through gg fusion, but this becomes a strong 
anticorrelation as the mass of the Higgs boson increases, as the gluons are in a similar x 
range as those responsible for ti production. Associated Higgs boson production {WhP) is 
strongly correlated with Z production for low Higgs masses but becomes decorrelated for 
higher masses.^ There is only a mild correlation between the production of a Higgs boson 
through vector boson fusion and Z production over a wide range of ruh- 



The scale dependence of the NLO tt cross section may be reduced by threshold resummation [72, 82, 83, 88], 
which includes higher-order logarithmic terms that enforce renormalization group invariance. In our study 
the scale dependence is viewed as an estimate of all NNLO contributions, including potentially sizable 
higher-order terms not associated with the threshold logarithms. A more conservative estimate of full 
NNLO effects is provided by the scale dependence of the fixed-order NLO cross section (11), rather than 
that of the threshold-resummed NLO cross section. 

Similarly, WhP associated production {mh < 200 GeV) and Z boson production are strongly correlated 
at the Tevatron. 
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Figure 15: The correlation cosine cos 99 for Higgs boson searches at the LHC with respect to Z boson 
production at the LHC (sohd) and Tevatron (dots), and tt production at the LHC (dashes), plotted 
as a function of Higgs mass. Separate markers denote correlations of W ^ t-channel single-top cross 
sections at the LHC and Z cross section at the Tevatron with respect to Z and tt cross sections at 
the LHC. 



The correlation curves with respect to it production basically form a mirror image to the 
previous curves, since ti pairs at the LHC are predominantly produced via gg fusion at large 
X. Thus, for example, the production of a Higgs boson through gg fusion goes from a mild 
anti-correlation with ti for low Higgs masses to a high correlation for large Higgs masses. 

PDF-induced correlations may follow a different pattern in other two-Higgs doublet mod- 
els (2HDM). For example, the MSSM Higgs boson production cs+c6 — > /i"*" is not particularly 
correlated with tt production at rrih ^ 500 GeV because of the uncorrelated contribution 
from the anti-strangeness PDF (cf. Fig. 15). The cs — > hf^ channel is absent in other 2HDM, 
such as the effective weak-scale 2HDM induced by top-color dynamics [120], where the scat- 
tering proceeds entirely via cb —>■ at the Born level. In those models the correlation of 
the ti cross section with the cross section is very strong at large nih (cosy? = 0.98 for 
ruh = 500 GeV). 

Higgs boson production at the LHC is affected by the new features of CTEQ6.6 PDFs in 
several ways. Table 5 lists the relative difference Aqm = o"6.i/o"6.6 ~ 1 between the CTEQ6.1 
and CTEQ6.6 cross sections, as well as CTEQ6.6 PDF uncertainties Apdf for Higgs boson 
production in a mass range rrih = 100 — 800 GeV. In processes dominated by light-quark 
scattering (vector boson fusion, W^h, Z^h, and Ah^), the most tangible differences between 
CTEQ6.6 and CTEQ6.1 (compared to the PDF uncertainty) occur at of order 100 GeV, 
where Agm is close in magnitude to Apcp, reflecting the enhancement in CTEQ6.6 u and d 
PDFs at X = 10~^— 10~^. In gluon-gluon fusion gg — > /i, the difference between CTEQ6.6 and 
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6.1 is well within the PDF uncertainty {Agm < Apof), although Aqm becomes comparable 
to ApDF for heavy Higgs masses {mh ^ 500 GeV). Similarly, Agm is smaller than Ap^p in 
heavy-quark scattering cb ^ h'^. The most striking differences occur in cs ^ h'^, because 
the CTEQ6.1 and CTEQ6.6 strangeness distributions disagree by a large amount in most 
of the X range. 

5. CONCLUSION 

The new CTEQ6.6 global analysis incorporates the latest improvements in the pertur- 
bative QCD treatment of s, c, and b quark PDFs introduced in Refs. [5-7]. It predicts 
substantial modifications (comparable to, or exceeding in magnitude the NNLO contribu- 
tions) in high-energy electroweak precision cross sections. Theoretical improvements of this 
kind must be accompanied by the development of efficient strategies to understand and 
reduce the remaining uncertainties in the PDF parameters. This work presents a novel cor- 
relation analysis, a technique based on the Hessian method that links the PDF uncertainty 
in a hadronic cross section to PDFs for physical parton flavors at well-defined {x, /i) values. 

We apply the correlation analysis to reveal and explain several regularities observed in 
the PDF dependence, such as strong sensitivity of W, Z production at the LHC and ti 
production at the Tevatron (processes dominated by qq scattering) to uncertainties in the 
gluon and heavy-quark PDFs; the leading role played by the strangeness distribution s{x, jj) 
in the PDF uncertainty of the ratio rzw = of the LHC Z and W cross sections; 

and intriguing PDF-induced anticorrelations between the LHC W cross sections and 
processes dominated by large- a; gluon-scattering, such as heavy Higgs boson production via 
gluon fusion. 

In our study we identify pairs of hadronic cross sections with strongly correlated or 
anticorrelated PDF dependence, i.e. with the correlation cosine cosip close to 1 or -1. 
Such pairs are especially helpful for constraining the PDF uncertainties, in view that a 
combination of the existing CTEQ6.6 and upcoming LHC constraints on one cross section 
in the pair is guaranteed to reduce substantially the PDF uncertainty on the second cross 
section [cf. the discussion accompanying Eq. (8)]. In addition, ratios of correlated (but 
not anticorrelated) cross sections have greatly reduced PDF uncertainty, as follows from 
Eq. (10). For this reason, it is beneficial to normalize an LHC cross section to a standard 
candle cross section with which it has a large PDF-induced correlation. 

We point out a potentially valuable role of precise measurements of tt cross sections at 
the Tevatron and LHC for constraining the gluon PDF at large x and normalizing the LHC 
cross sections that are anticorrelated with W and Z boson production. If both theoretical 
and experimental uncertainties on tt cross sections are reduced to a level of 3-5%, as may 
become possible in the near future, the tt cross sections would provide an additional standard 
candle observable with useful complementarity to Z and W boson cross sections. These 
measurements will be essential for reducing theoretical uncertainties in single-top and Higgs 
boson production, and for constructing cross section ratios with small PDF uncertainties. 
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Table 1: Total cross sections a, PDF-induced errors Aa, and correlation cosines cosip for Z , W , 
and tt production at the Tevatron Run-2 (Tev2) and LHC, computed with CTEQ6.6 PDFs. 
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Table 2: The fitting parameters (A,S,C) for the parametric form (11) of the CTEQ6.6M total 
cross section for inclusive tt and single-top production at the Tevatron and LHC, evaluated at NLO 
in the QCD coupling strength. Also shown are the relative scale and PDF errors, and I^pdf 
at mt = 171 GeV. 
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Table 3: Correlation cosines cos 93 between single-top, W, Z, and tt cross sections at the Tevatron 
Run-2 (Tev2) and LHC, computed with CTEQ6.6 PDFs. 
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Table 4: CTEQ6.6M total cross sections a and PDF errors Act for Higgs boson production at 



= 14 TeV shown in Fig. 15. 
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Table 5: Relative differences Aqm = o'e.i/o'e.e — 1 between CTEQ 6.1 and CTEQ 6.6 cross sec- 
tions for Higgs boson production at the LHC, compared to the PDF uncertainties Apop in these 
processes. The Ah^ cross section is for combined production of positively and negatively charged 
Higgs bosons, with m/j being the mass of the CP-odd boson (m^ = itia), and m/j± given by 
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